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Abstract 

This article is devoted to developing a theory for effective kernel interpolation and 
approximation in a general setting. For a wide class of compact, connected Rieman- 
nian manifolds, including the important cases of spheres and SO{3), we establish, using 
techniques involving differential geometry and Lie groups, that the kernels obtained as 
fundamental solutions of certain partial differential operators generate Lagrange functions 
that are uniformly bounded and decay away from their center at an algebraic rate, and 
in certain cases, an exponential rate. An immediate corollary is that the corresponding 
Lebesgue constants for interpolation as well as for L2 minimization are uniformly bounded 
with a constant whose only dependence on the set of data sites is reflected in the mesh 
ratio, which measures the uniformity of the data. The kernels considered here include the 
restricted surface splines on spheres, as well as surface splines for 5*0(3), both of which 
have elementary closed-form representations that are computationally implementable. In 
addition to obtaining bounded Lebesgue constants in this setting, we also establish a "ze- 
ros lemma" for domains on compact Riemannian manifolds - one that holds in as much 
generality as the corresponding Euclidean zeros lemma (on Lipschitz domains satisfying 
interior cone conditions) with constants that clearly demonstrate the influence of the ge- 
ometry of the boundary (via cone parameters) as well as that of the Riemannian metric. 



1 Introduction 

Radial basis functions (RBFs) have proven to be a powerful tool for analyzing scattered data 
on i?". More recently, spherical basis functions (SBFs), which are analogs of RBFs on the 
n-sphere, and periodic basis functions (PBFs), which are analogs of RBFs on the n-torus, 
have had comparable success for analyzing scattered data on these manifolds. A theoretical 
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drawback is that most RBFs are globally defined, thin plates splines are an example, and even 
those that are locally defined such as Wendland functions behave globally when approximating 
at densely packed data sites. Nevertheless, certain RBF approximants, in their numerical 
implementation, exhibit localized behavior, i.e., changing data locally only significantly alters 
the interpolant locally. Since the pioneering work of Duchon O [8], it has been a mystery why 
RBF/SBF approximants display local behavior even though the bases are globally supported. 
It was long suspected that there were "local bases" hidden within the space of translates of 
RBFs/SBFs. 

A major objective of this paper is to establish that, for spheres and S0{3), there are 
closed-form kernels whose associated approximation spaces possess highly localized bases, in 
the form of Lagrange functions for given scattered data. Our previous work [16j established 
the existence of such bases on compact manifolds, but the kernels we constructed did not have 
closed form. 

We will carry out the construction and proofs that establish the existence and properties of 
these closed- form kernels in the context of of more general manifolds, thus achieving another 
objective: obtaining results for a broader class of kernels on manifolds than the ones we 
treated in [16], where only reproducing kernels for Sobolev spaces were dealt with. 

We also address similar issues for conditionally positive kernels on manifolds, where a 
given space of functions is to be reproduced. In the case of M*^, this involves little more than 
polynomial reproduction. For manifolds - even for spheres and 5*0 (3) - the kernels and spaces 
are not so simple, and new techniques are required to deal with the problem. 

Goals Given a manifold M, a finite set of points X = {xi, . . . ,xn} C M and a kernel 
k: M X M — )• M, one may attempt to fit a smooth function using functions from Vx ■= 
span{fc(-, Xj), Xj € X}, or more generally, to use functions of the form 



where the supplementary function p comes from a simple space (like polynomials or spherical 
harmonics). The framework described above applies to fitting data by means of interpolation, 
least squares or near interpolation with a smoothing term. 

This article is devoted to developing a theory for effective kernel approximation in a general 
setting. The problem is described as follows: we seek kernels k: M x M — t- M for which interpo- 
lation is well posed and that have a convenient closed form representation allowing for effective 
computations. Furthermore, we are interested in aspects of the interpolants/approximants 
concerning stability, locality and so on. 

In [16j and |l5j, we developed a theory for compact Riemannian manifolds using positive 
definite "Sobolev kernels." The theory developed there addresses and answers questions con- 
cerning properties of bases for Vx, properties related to locality, stability of approximation 
and interpolation, and other matters. In this paper, these questions, which are listed below, 
are addressed and answered for a broad class of kernels on M that are Green's functions for 
certain elliptic operators, and, when the manifold is a sphere, real projective space or SO{3), 
are computationally implementable as well. 



N 





2 



Locality. Are there local bases for Vx? That is, are there bases similar to those for 
wavelet systems or B-splines [Gj Chapter 5]? At a minimum, we would like a basis {vj} to 
satisfy < r{d{x, Xj)), with r a rapidly decaying function. 

Lp conditioning. Are there bases that are well conditioned in Lp, after renormal- 
ization? That is, can we find bases for which there constants ci,C2 such that ci||a||£p < 

II Z^jLi OjfjIiLp < C2||a||^p, with ci,C2 independent of A^, and, after a suitable normalization, 
independent of p7 

Marcinkiewicz-Zygmund property. Does the space Vx possess a Marcinkiewicz- 
Zygmund property relating samples to the size of the function? For s £ Vx, this means 
that the norms \\j i— )• s(a;j)||£p and ||s||ip are equivalent, with constants involved independent 
of AT. 

Stability of interpolation. Is interpolation stable? Is the Lebesgue constant bounded 
or, more generally, is the p norm of the interpolant controlled by the £p norm of the data? 

Stability of approximation in Lp. Is approximation by L2 projection stable in Lp? 
Here, 1 < p < 00. In particular, what we want is that the orthogonal projector with range Vx 
be continuously extended to each Lp, and that it has bounded operator norm independent of 
A^. 

The Sobolev kernels we dealt with in [16] do not possess simple, closed form representa- 
tions, even when the underlying manifold is a sphere; they are defined indirectly, as reproduc- 
ing kernels for certain Sobolev spaces. To applied effectively to data fitting problems, such a 
kernel should have an implementable characterization, by which interpolation, approximation 
or other computational problems can be treated. In the important cases relating to spheres 
and S0{3), we exhibit computationally implementable kernels. In particular, these kernels 
include restricted surface splines on spheres, and surface splines on SO (3), both of which 
have simple closed form representations. Furthermore, for both of these cases, theoretical ap- 
proximation results concerning direct theorems, inverse theorems, and Bernstein inequalities 
are known to hold [22]. In conjunction with the stability of interpolation and least-squares 
approximation in Lp, both yield new, precise error estimates for these implementable schemes. 

Kernels The class of kernels considered in this paper are those kernels k that act as fun- 
damental solutions for elliptic differential operators of the form Cm = 11 ~ '"i) lower 
order perturbations of these. The origin of this approach lies in the work of Duchon, [9l [8] on 
surface splines, where the underlying kernel is the Green's functions for iterated Laplacian, 
A™, on M*^. Such kernels have also been used on Riemannian manifolds, [101 127] . For this 
reason, we call them polyharmonic; see Definition 13.21 Throughout this article, we use hfYi to 
denote a generic polyharmonic kernel. 

This is a classic family of kernels and is sufficiently robust to include many interesting 
examples. For instance, such kernels have also been in use for some time on spheres, and have 
formed one of the earliest families of SBFs (see fTT] and references for a list of early examples) . 
In this setting, certain careful choices of these kernels result in the complete family of surface 
splines restricted to the spher^, introduced in [24] , which we define below in (13. 2p and denote 
in the "zonal" form as km{x, t) = (j)s{x ■ t). Here s is related to m via m = s + d/2. 

related problem, which can be considered a generalization of this particular set up has recently been 
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It also includes the surface splines on iS'0(3), introduced in [T7] and defined below in (j3.4p . 
and denoted throughout the paper by k^- 

A second type of kernel, ostensibly different from the polyharmonic kernels, are the Sobolev 
(or Matern) kernels, which have been introduced for compact Riemannian manifolds in |16] . 
These kernels come about as reproducing kernels for Sobolev spaces. We denote such kernels 
by Km{x,y), where m indicates the order of the Sobolev space. A corollary of the results 
presented in this paper is that, in many cases, the Sobolev kernels are in fact polyharmonic 
kernels. There is an operator Cm for which Km{x,y) is the fundamental solution. 



Kernel 


Notation 


Location in manuscript 


Restricted surface spline 


{x,y) ^ (psix- y) 


Example 2 in 13.21 


Surface spline on S0(3) 




Example 3 in 13.21 


Polyharmonic kernel 


k 


Definition [32] 


Sobolev kernel 




tl6i 3.3] 



Table 1: Index for kernels used 



Outline The layout of this paper is as follows. Following the introduction and background. 
Section [2] deals with certain geometric notions relevant to this article. Section [3] treats inter- 
polation by conditionally positive definite kernels, the function spaces that they generate, and 
the nature of their interpolants. We discuss some important examples on well known mani- 
folds, including spheres and the rotation group. Finally we define precisely the polyharmonic 
kernels, which are the kernels that we treat in our main results; they include the examples 
we provided earlier. We demonstrate that they are conditionally positive definite, identify 
the seminorm of the native spaces associated with these kernels, and discuss the variational 
problem associated with their interpolants. 

The relationship between the polyharmonic kernels and the Sobolev kernels of [16j will be 
covered in Section We show that, under certain conditions, the polyharmonic operators 
Cm can be expressed as combinations of operators generated by covariant derivatives, and 
vice versa (this is done in Lemma 14. 3p . This allows us to conclude that Sobolev kernels 
are examples of polyharmonic kernels. On the other hand, it permits us to demonstrate, in 
Section 14.21 that the native space seminorms associated with polyharmonic kernels exhibit 
the same behavior (metric equivalence to Euclidean Sobolev seminorms, zeros lemmas, etc.) 
as the native space norms associated with Sobolev kernels. 

The main results of the paper are given in Section [5j Namely, the Lagrange function 
associated with a kernel km is rapidly decaying, and gives rise to a uniformly bounded Lebesgue 
constant and a uniformly bounded L2 minimization projector. The properties mentioned 

considered by Fuselier and Wright, [12]. There, kernel interpolation is considered on manifolds that are embed- 
ded in R"* by using the restriction of various other RBFs to the manifold - this is accomplished by constructing 
interpolants in the ambient Euclidean space and then restricting these to the manifold. (In contrast, we work 
directly with the manifold, making use of its intrinsic structure.) 
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above - concerning locality, stability, conditioning, and so on - then follow immediately. We 
then discuss implications of this for surface spline kernels on spheres and SO (3). 

Essential to our proofs in Section [5] are theorems giving Lp Sobolev-space estimates for 
functions having zeros quasi-uniformly distributed on a domain fi, with 90 being Lipschitz. 
Such theorems may hold both in and on M itself, and in Section A, we treat both cases. For 
the case of a manifold M, these theorems involve geometric ideas; in particular, they require 
use of a minimal e-set of points in M (cf. [13j), which replaces a simpler set in W^. The results 
for M turn out to be intrinsic and hold in the same generality as those in the Euclidean case. 
The bounds and the condition on the meshnorm reflect geometric properties - particularly, 
parameters from the cone condition on dQ, properties of the manifold M, and parameters of 
the Sobolev spaces themselves - but are independent of the volume and diameter of Q. 

2 Background 

We now discuss some relevant details about analysis on compact, complete, connected C°° 
Riemannian manifolds. This is the same setting as [16j. Refer to it for a more detailed 
treatment and further references. 

Throughout our discussion, we will assume that (M,^) is a d-dimensional, connected, 
complete C°° Riemannian manifold without boundary; the Riemannian metric for M is g, 
which defines an inner product gp{-, •) = (•, ■)g^p on each tangent space TpM; the corresponding 
norm is | -{g^p- As usual, a chart is a pair (il, </>) such that il C M is open and the map : il — )• 
is a one-to-one homeomorphism. An atlas is a collection of charts {{ila,4>a)} indexed by a 
such that M = Ua^a and, when cpai^a) H (/)/3(il^) 7^ 0, (/>^ o is C°°. In a fixed chart (il, (p), 
the points p E il are parametrized by p = </)~^(x), where x = {x^, . . . , x*^) E C/ = </)(il). 

In these local coordinates, TMp and T*Mp, the tangent and cotangent spaces at p, have 
bases comprising the vectors ej = (gfj)^, j = 1 ■ ■ - d and e'^ = {dx^)p, k = 1 . . . d, respectively. 

These vary smoothly over U = (/>(il) and form dual bases in the sense that e'^(ej) = = 6j. 
In the usual way, the inner product (•, ■)g^p provides an isomorphism between the cotangent 
and tangent spaces. Thus, regarding e*^'s as vectors, we have that {e^,ej)g^p = 5^. A vector 
V can be represented either as v = "^jV^ej or as v = X^fc^fc^'^i ^"''^ and VkS are the 
contravariant and covariant components of v, respectively. Relative to these bases, the inner 
product (•, ■)g^p has the form 

d d 

{u,\r)g^p = ^ giju'v^ = ^ g'^UiVj , where gij = {ei,ej)g^p and g'^ = {e\e^)g^p (2.1) 

i,j=l i,j=l 

The matrices (gij) and {g^^) are inverse to each other, and are of course symmetric and positive 
definite. The inner product (v, w)g^p is itself independent of coordinates. In addition, if v 
and w are C°° vectors fields in p, then {v,w)g^p is also C°°. 

The metric g also induces an invariant volume measure d/x on M. The local form of the 
measure is d^{x) = \J A.Q\,{^g)dx^ ■ ■ ■ dx*^, where det{g) = det{gij). 

Geodesies are curves 7 : M — )• M that locally minimize the arc length functional, \'y\g^pdt. 
If we use the arc length s as the parameter (i.e., t — )• s), then, in local coordinates, a geodesic 
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satisfies the Euler-Lagrange equations: 



Vr^^— -0 where r^-- V , ^ _ 2) 

«J=1 me{l...d} 

The are the Christoffel symbols. 

A geodesic solving (12. 2p is specified by giving an initial point p G M, whose coordinates we 
may take to be (x^(0), . . . ,x'^(0)) = 0, together with a tangent vector tp having components 
^(0). A Riemannian manifold is said to be complete if the geodesies are defined for all values 
of the parameter s. All compact Riemannian manifolds without boundary are complete, and 
so are many non-compact ones, including R'^. 

We define the exponential map Exp^ : TpM — )• M by letting Expp(O) = p and Expp(stp) = 
7p(s), where 7p(s) is the unique geodesic that passes through p for s = and has a tangent 
vector 7p(0) = tp of length 1; i.e., (tp,tp)gp. = 1. 

Although geodesies having different initial, non-parallel unit tangent vectors tp = 7p(0) 
may eventually intersect, there will always be a neighborhood itp of p where they do not. In ilp, 
the initial direction tp, together with the arc length s, uniquely specify a point q via q = 7p(s), 
and the exponential map Expp is a diffeomorphism between the corresponding neighborhoods 
of in TpM and p in M. In particular, there will be a largest ball -6(0, rp) G TpM about the 
origin in TpM such that Expp : B{0, rp) — )■ b(p, rp) C M is injective and thus a diffeomorphism; 
rp is called the injectivity radius for p. By choosing cartesian coordinates on 5(0, rp), with 
origin 0, and using the exponential map, we can parametrize M in a neighborhood of p via 
q = Expp(x), X G TpM. 

The injectivity radius of M is rM := iufp^Mfp- If < tm < oo, the manifold is said 
to have positive injectivity radius. For any r < tm and any p G M, the exponential map 
Expp : B{0,r) — t- h{p,r) is a diffeomorphism. 

We make special note of the fact that, for a compact Riemannian manifold, the family of 
exponential maps are uniformly isomorphic; i.e., there are constants < Fi < 1 < r2 < oo so 
that for every pQ £ M and every x,y £ i?(0,r), where r < tm/S, 

ri|x-y| < dist(Exppj,(x),Expp^(y)) <r2|x-y|. (2.3) 

When we use Fi and F2 in this paper, we always assume that there is some fixed radius 
smaller that r^i on which they are computed. This avoids a tiresome repetition of this fact 
throughout the paper. 

An order k covariant tensor T is a real-valued, multilinear function of the /c-fold tensor 
product of TpM. We denote by TpM the covariant tensors of of order k at p. In terms of 
the local coordinates, there is a smoothly varying basis e*i iS" • • • (8) e**= for the A;-fold tensor 
product of tangent spaces. Thus, the covariant tensor field T of order A; on il can be written 
as 

ie{i...d}'' 
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where we adopt the convention i = (zi, . . . The Tj are the covariant components of T. 
The metric gij is itself an order 2 covariant tensor field. One can also define contravariant 
tensors and tensors of mixed type. 

Because T^M = TpM (8) • • • (8> TpM (k times), the metric g induces a natural, useful, 
invariant inner product on T^M; in terms of covariant components, it is given by 

(S,T)3,p= g''''---9''"'SiT^. (2.4) 

i,ie{i...ci}* 

The covariant derivative^ or connection^ V associated with (M,(7) is defined as follows. If 
T is an order k (covariant) tensor, then the covariant derivative of T is 

je{i...d}ie{i...dY r=i se{i...d} ^ 

(VT)s,j 

which is an order k + \ covariant tensor with components (VT)jj. The T^j are the Christoffel 
symbols defined earlier. 

A smooth function / : M — )• M is a order tensor and so V/, which is the gradient of / is 
an order 1 tensor, V^/ is an order 2 tensor. Continuing in this way, we may form V'^/, which 
is an invariant version of the ordinary A;*^ gradient of a function on M*^. (The superscript k 
here is an operator power, not a contravariant index.) The components of the k*^ covariant 
derivative of / have the form 

(VV(x)).- = (aV(x)).' + E E (2-5) 

m=l je{l...d]^- 

where 

and where the coefficients x i— )• Al{x) depend on the Christoffel symbols and their derivatives 
to order A; — 1, and, hence, are smooth in x. This can also be written in standard multi- 
index notation. Let ai, 0:2, . . . ,ad be the number of repetitions of 1, 2, . . . , d in j, and let 
a = (ai , . . . , ad) ■ Then, 

^^^^^-^ ■= d{x^r^ • • • dix'^r^ ^ ° =■ ° '^''^ H = E«^ = -- M 

Another important quantity that we need to deal with is the adjoint of the covariant 
derivative V*. This operator is defined by Jj^(V*T, S)g^p(i^ = f^{T,VS)g^pdn, where the 
inner product is given by (|2.4|) . and it takes an order k + 1 tensor to an order k tensor. The 
coordinate form of V*T is obtained via integration by parts: 

(V*T)5 = -^5^-^(VT)5,,-fc (2.7) 

id 



7 



We can combine covariant derivatives and their adjoints to get scalar operators. In par- 
ticular, if / : M ^ M is C°°, then V^f is an order k tensor. By applying V* to it, we get 
(V*)''V*^/, which is a scalar. (Note that (V*)'' = (V'')*.) The Laplace-Beltrami operator 
A := — V*V. In coordinates, again letting det{g) = det{gij), we have that 



Sobolev spaces on subsets of M Sobolev spaces on subsets of a Riemannian manifold 
can be defined in an invariant way, using covariant derivatives [I]. In defining them, we will 
need to make use of the spaces Lp, Lq. To avoid problems with notation, we will use the 
sans-serif letters p, q, rather than p, q, as subscripts. Here is the definition. 

Definition 2.1 ([1, p. 32]). Let C M 6e a measurable subset. For all 1 < p < oo, we define 
the Sobolev space W^{n) to be all f : M R such that, for < k < m, \'^^f\g,p in Lp(fi). 
The associated norms are as follows: 

_f {ET=oIn\^'f\lpMp))'^\ l<p<oo; 



maxo<A;<m|||V /|g,p||i^(Q), P 



oo. 



When p = 2, the norm comes from the Sobolev inner product 

m 

{f,9)m,n ■■= {f,9)wrin) - E / (^V, V^^) dfi{p). (2.9) 



A;=0' 



We also write the p = 2 Sobolev norm as ||/||^ q (/' f)m,n- When Q = Wl, we may suppress 
the domain: {f,g)m = {f,g)mM and 



Metric equivalence The exponential map allows us to compare the Sobolev norms we've 
just introduced, to standard Euclidean Sobolev norms as follows: 

Lemma 2.2 Lemma 3.2]). For m G N and < r < rM/3, there are constants < ci < C2 
so that for any measurable Q C B^, for all j € N, j < m, and for any pQ G M, the equivalence 

holds for all u : 'Expp^(il.) — t- M. The constants ci and C2 depend on r, m and p, but they are 
independent of and p^ . 



Besov spaces on M Besov spaces can be defined and characterized in many equivalent ways. 
For a discussion, see Triebel's book |281 1.11, and Chapter 7] and the references therein. Our 
definition follows Triebel's. 
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Definition 2.3. For < s < m and 1 < p < oo, we define the Besov space ^^^^(m) as the 
collection of functions in Lp(M) for which the following expression 

\\f\\Bs^(M) ■.= snpt-^K{f,t) 

t>0 

is finite, where the K -functional K(f, •) : (0, oo) — )• (0, oo) is defined as 

K{f, t) := inf { 11/ - gU, + t'^'"||5||^™(M) : 9 ^ W^{M)} . 

For p = oo, the definition is the same after substituting Lp(M) by C(M) and VF™(M) by 
C'"(M). 

2.1 Notation 

In order to distinguish balls on M'^ from tliose in M, we denote the ball centered at p G M 
having radius r by h{p,r). (Euclidean balls are denoted B(x,r).) Given a finite set S C M, 
we define its mesh norm (or fill distance) h and the separation distance q to be: 

/i := sup dist(p, S) and q := inf dist(^, C)- (2-10) 

The mesh norm measures the density of H in M, the separation radius determines the spacing 
of H. The mesh ratio p := h/q measures the uniformity of the distribution of H in M. If p is 
bounded, then we say that the point set H is quasi-uniformly distributed, or simply that H is 
quasi-uniform. 

3 Interpolation by Kernels 

The purpose of this section is to discuss further this interpolation problem and to present 
the kernels we employ. The kernels we consider are fundamental solutions of certain elliptic 
PDEs. They happen also to be conditionally positive definite, a well known class for which 
interpolation is understood. In particular, interpolation is well posed, and has a dual nature, 
as best interpolation from a function space. 

In 13. II we discuss interpolation with conditionally positive definite kernels, and present the 
associated problem of best interpolation. In 13.21 we present some motivating examples for our 
problem: surface spline interpolation on spheres and on S0{2>). In 13.31 we give the formal 
definition of the kernels we use and the operators they invert; we also discuss the associated 
variational problem they solve. 

3.1 Interpolation with conditionally positive definite kernels 

The kernels we consider in this article are conditionally positive definite on the compact 
Riemannian manifold. As a reference on this topic, we suggest |101 Section 4]. 
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Definition 3.1. A kernel is conditionally positive definite with respect to a finite dimensional 
space n if, for any set of centers H, the matrix [k{^, ()^^ is positive definite on the subspace 

of all vectors a G satisfying Xl^es '^(.PiO — f'^^ p Gil. 

This is a very general definition which we will make concrete in the next subsections. 
Given a complete orthonormal basis (^j)jgN) of continuous functions, normalized in (i.e., 
||(/)j||oo = 1) any kernel 

with coefficients k G ^i(N) for which all but finitely many coefficients k{j) are positive is 
conditionally positive definite with respect to IIj- = span(^j \ j E J), where J = {j \ k{j) < 
0}, since, evidently. 



provided ^^Oi^4'jiO = foi" J satisfying k{j) < 0. 

In this case if the set of centers H C M is unisolvent with respect to XIj- = span(99j | j G J') 
(meaning that p ElLj and p{^) = for ^ G S implies that p = 0) then the system of equations 



E^es«e<^j(0 = 



has a unique solution in C" x C"^ for each data sequence (y^)^^^ G C". When data is sampled 
from a continuous function at the points S (i.e., j/^ = /(C)) this solution generates a continuous 
interpolant: 

isf = ik,j,Ef = E «?^(-' + E ^j^j' 

with the property that it is the minimizer of the semi-norm ||H||j. j-, called the "native space" 
norm, given by 



k,J 



E 



ml 



(3.1) 



over all functions u = Ylijen'^ij) which n(^) = y^, ^ G H. If k is conditionally positive 
definite with respect to the set 11 j-, it will be conditionally positive definite with respect to 
n ji for any J' D J. For this reason, the interpolant and norm both are decorated by k and 
J. 

This has the consequence that any two conditionally positive definite kernels fc, k' which 
have eigenfunction expansions that coincide on all but finitely many indices (say I), produce 
the same interpolants. That is: /fc,x,s = -ffe',x,s- 
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3.2 Examples of conditionally positive definite kernels 

Example 1 (Surface Splines) . As a first example of a conditionally positive definite kernel, we 
take M = M"^, and consider the kernels km{x, a) = (f)g{x — a) given by the functions 



= C, 



rriA 



\x\^^ log |x| dis even 
IxP* dis odd 



where s + d/2 = m. For a certain Cm,d 7^ 0, this is a fundamental solution for the operator 
A™. 

Because of the positivity of the generalized Fourier transform, one can see that (f)g is con- 
ditionally positive definite on with respect to Hm-i- These have been studied by Duchon 
[9l |8], and they comprise some of the earliest and most popular examples of conditionally 
positive definite kernels. 

Although our focus in this paper is on kernels on compact manifolds, the family of surface 
splines acts as a useful benchmark, since they have a simple, direct representation, as well as 
being conditionally positive definite, not to mention that for certain interpolation problems, 
their Lagrange functions decay rapidly (this was demonstrated in a least squares sense by 
Matveev in [211 Lemma 5] and pointwise in [1^) and have a uniformly bounded Lebesgue 
constant (cf. [16]). 



Example 2 (Restricted Surface Splines on S'^). When M = S'^, the eigenvalues of the Laplace- 

:d+e)r{e+d- 
r{<?+i)r(d) 



Beltrami operator are fj,i = i{£ + d — 1) and each eigenvalue has N{d, i) — (^^+^)^(^+'^ ^) 



linearly independent eigenfunctions, the spherical harmonics l^.m- 

We now introduce a family of kernels known as the restricted surface splines. These are 
kernels indexed by m G N, m > d/2. By writing m = s + d/2, we give the zonal expression 

^.(t):=|"^-*"^^°^"^-*" (3.2) 

When d is even, s is integral and the first formula is used. When d is odd, the second is 
used. For a given d and an integer m > d/2, we write km{x,y) = (j)s{x ■ y) to denote the 
corresponding kernel on S"^. 

A spherical harmonic expansion of the restricted surface splines can be found in ^ Equa- 
tions (2.12) & (2.20)]. It is known that, for m > d/2, kmix, y) = T,j km{£, j)yi,j{x)Yij{y), 
where 



+ (^'lV-(--l/2)= 



for i> s. (3.3) 



kmi^yj) — Cm 

When d is odd, this equation holds for all i. 

From this formula, it follows that km is conditionally positive definite with respect to the 
space Ili^gj := span(y£j \ £ < s, j < N{d,£)). 

A second consequence is that, by a possible slight correction of the spherical harmonic 
expansion (discussed below), km is the fundamental solution for a differential operator of 
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order 2m that is polynomial in A: 

m 

J^m ■■= CmlliA-iu- d/2){v + d/2 + 1)] . 

We note that when d is odd, the operator Cm is invertible on W2^{E>'^). Indeed, it is 
nonvanishing on each spherical harmonic Yi ^- 

When d is even, the Fourier coefficients of the kernel follow (j3.3p for £ ^ s only, but 
annihilates spherical harmonics of degree s or less. Indeed, in this case, we can re-index the 
operator to get: 

d/2-l m 

Cm = Cm H [A - - d/2){u + d/2 - 1)] J] [A - - d/2){u + d/2 - 1)] 

v=l v=d/2 
d/2-l m-d/2 

= CmW [A-{v-d/2){u + d/2-l)] W [A- J(J + d-l)]. 

u=l J=0 

So Cm annihilates all the spherical harmonics of order up to s = m — d/2. 

In other words, for sufficiently smooth functions, say / G C^™(§'^) represented by the 



f{x)= Cm[f]{a)(l)s{x ■ a)dfi{a) +pf 



where we add a spherical harmonic term pf = J2t=o'^f=i'^^ ^ when d is even 

(when d is odd, pf = 0). 

Example 3 (Surface Splines on SO (3)). When M = 50(3), the group of proper rotations of 
M^, the eigenvalues of the Laplace-Beltrami operator are /u^ = £{£ + 1) and each eigenvalue 
is associated with A^(^) = (1 + 2£)^ linearly independent eigenfunctions, called Wigner D- 
functions and denoted by (-^^j J(|j|,|;,|<£)- 

For m >2 and s = m — 3/2, the surface spline kernels are 

k^(x,,)=fsmf^^^lV'"", (3.4) 



where uj{z) is the rotational angle of z G S0{3). The corresponding Wigner D-function 
expansion km(x,y) = Yle'l^j i,^rni£,j, '-)Dj^}ix)Dj^j{y) is discussed in [17, Lemma 2], where 
it is shown that for some Cm 7^ 0, 



kmii,j,i) =Cm 

i/=— (m— 1) 



+ ^+2 



Thus, km is conditionally positive definite with respect to the space 11^-2 = span{Dj^^ \ £ < 
m-2,\j\,\i\<£}. 



12 



It also follows (from [17^ Lemma 3]) that is the fundamental solution for the differential 
operator of order 2m having the form: 



m—l 



C^ll[A- {u' - 1/4)] 



in the sense that for f €z C' 



fix) 



f 




holds true. 

3.3 Poly harmonic and related kernels 

The kernels we wish to treat are fundamental solutions of differential operators that are 
polynomial in the Laplace-Beltrami operator, or are directly related to them. Since, on 
a compact Riemannian manifold A is a self adjoint operator with a countable sequence of 
nonnegative eigenvalues Xj < Aj+i having +oo as the only accumulation point, we can express 
the kernel in terms of the associated eigenfunctions A(f)j = \j4>j- We now make this clear 
with a formal definition. 

Definition 3.2. Let m G N such that m > d/2. We say that the kernel ■ M x M ^ M. is 
polyharmonic if the following hold: 

1. There exists a polynomial Q{x) = Y^^= Hfni^), with the highest order coefficient 

Cm > 0, so that Q{x) > for all x sufficiently large. Let the corresponding differential 
operator of order 2m > d be given by 



and let J dl^ be a finite set that includes all j for which the eigenvalue Q{Xj) of Cm 
satisfies Q{)^j) < 0. (In addition to this finite set, J may also include a finite number 
j 's for which Q{Xj) > 0.) 

2. The kernel has the eigenfunction expansion km{x,y) = J2jeN^m{j)ipj{x)ipj{y), with 
coefficients km{j) = l/Q{Xj) for j ^ J'. (On J , kmij) can assume arbitrary values.) 

It follows immediately from this definition that km is conditionally positive definite with 
respect to the finite dimensional space 11 j- = span^gj- (^j. Another consequence is that, for 



m 



/:^ = ^c,A^ = Q(A) 



'OO 




(3.5) 



where pj = Sjgj-projj/ is the orthogonal projection onto IIj-. 
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As previously stated, the interpolation operator Is,km,J produces the minimizer of the 
seminorm |||u|||^^ Since km{x,y) = Y.jijkm{j)Vj{x)^j{y) and, for j ^ J, 



\ -1 



which is the inverse symbol of Cm, it follows from (j3.ip that 

3tJ ^^^^ j€J 



k{j) ^ 



This relation connects the norm j with the quadratic form {CmU, u)^2(M)- The goal of 

the next section is to study this quadratic form. 



4 Operators and quadratic forms 

Of the two quadratic forms considered, the one derived from the native space seminorm: 
lll^lllfcm derived from the operator [u]^ := {CmU, u)2,2(m)) the latter has much to 

offer from the point of view of analysis, but the former is tied to the variational problem sat- 
isfied by the kernel interpolants. The object of this section is to attain a better understanding 

of 

To this end, we seek an analogue of the bilinear form (£n, v)^^^^-^ - one that is defined on 
measurable subsets of M. A reasonable goal would be to find a form that is comparable to 
the corresponding Sobolev form YlY=o('^''^h,^^ where 

{u,v)j^Q:= / {V^u,V^v)g^r,dij{x). 
Jn 

This is the bilinear form used to define the Sobolev space inner product: (I2.8p of Definition 
Ofor ncM. 

The rest of this section is structured as follows. In l4.1l we demonstrate that on a wide class 
of manifolds, the elliptic operator composed of covariant and contravariant derivatives, which 
is at the heart of [16] . is a polynomial in A and, conversely, the Laplace-Beltrami operator 
has an expansion in terms of these elliptic operators. This permits us immediately to classify 
the Sobolev kernels on spheres (as well-known kernels of a type studied in [141 122j ) and to 
give concrete approximation results for them. In 14. 21 we present analogues to the bilinear form 
generated by Cm on measurable subsets. Using this, we demonstrate that this bilinear form 
behaves like a norm for functions with many zeros. 



4.1 The Laplace-Beltrami operator and the covariant derivative 

Simply considering {Cu,v) for measurable subsets C M is not suitable, since there will 
be many functions for which Cu may vanish on Q,. This is true even on M*^ when C = A*". 
In this case there are many polyharmonic functions (and even harmonic functions!) on a a 
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given subdomain Q that may have nonzero Sobolev norms, despite the fact that they are in 
the kernel of C. 

Guided by the observation that on M^, f^a vA'^u = J^di^^^vV^u) holds for test functions 
u,v, we first attempt to compare A*", the principle part of >C, to (V™)*V'". It is important to 
stress that (V™")* means the adjoint to V™, in the L2(M) inner product, as defined in (j2.7p . 
H To this end, we make the following assumption 

Assumption 4.1. For all A; € N, there exists a real polynomial Pk~i of degree k — 1, such 
that 

A class of Riemannian manifolds that satisfy this are the two-point homogeneous spaces 
[18| . both compact and non-compact. A manifold M is homogeneous if M = G/K, where G 
is a Lie group and K is a Lie subgroup of G. Two-point homogeneous means that for any two 
pairs of points p,q and p',q' such that the distances d{p,q)= d{p',q') there is an isometry 
$ G G such that p' = ^{p) and q' = ^{q). These manifoldqj have been completely classified 
(see \18i, p. 167 Sz p. 177] for lists), and include S'^ and the real projective spaces P"^. (The 
rotation group SO {3) = P^.) 

Lemma 4.2. Let Wl be a two-point homogeneous space. Then M satisfies A s sumption \4 ■ 1\ 

Proof. The proof proceeds in two steps. The first is showing that if <I> : M — )• M is a diffeo- 
morphism that is also an isometry (i.e., preserves distances), then the operator D := (V'^)*V^ 
is invariant in the sense that for any smooth function / : M — )■ M, Df = {D{f o <^)) o $^1. We 
will follow a technique used in [18, Proposition 2.4, p. 246]. Let (il, (/>) be a local chart, with 
coordinates x^ = (l>^{p), j = 1, . . . ,n for p G il. Since <I> is a diffeomorphism, (<5(il),<^ o <I)~^) 
is also a local chart. Let ip = cpo and use the coordinates = ip^{q) for q £ '^(H). The 
choice of coordinates has the effect of assigning the same point in M" to p and q, provided 
q = ^{p) - i.e., x^{p) = y^{q). Thus, relative to these coordinates the map <I> is the identity, 
and consequently, the two tangent vectors {-^)q G TqWi and {^j)p G TpM are related via 

So far, we have only used the fact that $ is a diffeomorphism. The map ^> being in addition 
an isometry then implies that 

\dy^ ' dy^ / ^^^^ \ ^ \dx^ J ' ^ \dx^ ) / q,{^p) \dxi' dx^ / ^' 

The expression on the left is the metric tensor at ^{p), g-'^{y), and on the right, the metric 
tensor at g^^{x). The equation above implies that, as functions of y and x, g^^{y) = g^^{x). 

^and not with respect to 1.2(0) - e.g., even though A™ = (-1)'"(V'")* V™ holds on R'', it is not the case 
that (-1)'" vA"^u = V^wV^-M for subsets n C R''. 

■^We note that they have also appeared in other approximation theory hterature. See, e.g., [3] 120) 
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This means that the expressions for the Christoffel symbols, covariant derivatives and various 
expressions formed from them wih, as functions, be the same. Since the operator D = {S/^)*\/^ 
is constructed from such objects, it fohows that Df = (D(/o<I)))o<i)~^, and so D is invariant. 

The second step makes use of two-point homogeneity. Since (V'^')*V'^ is invariant under 
every isometry <I> in G, applying [181 Proposition 4.11, p. 288] yields the result that (V'^)*V'^ 
is a polynomial in the Laplace-Beltrami operator: (V'^)*V'^ = a^+i^^ + afcA'^"^ + •••ao. 
Comparing terms in the highest order derivatives involved in coordinate expressions for both 
sides shows that a^+i = (—1)'^. □ 

Induction ensures that 
From this we have the following. 



Lemma 4.3. Suppose M. is a Riemannian manifold satisfying Assumption \4-l\ If Q{x) = 
CmX'^ + • • • + Co is a (real) polynomial of degree m, then there exist real numbers aj, with 

O-m = Cm, so that 

m 
j=0 

Conversely, for any constants bj, there is a real polynomial p for which the operators p(A) 
and '^j^Qn bj(V^)*V^ coincide. 

An immediate consequence is that the Sobolev kernels k^^m considered in [16j and |15j 
are Green's functions for operators of the form (5(A), with Q a real polynomial of degree m. 

We note, furthermore, that because the lead coefficient Cm of Q is assumed positive (see 
Definition 13. ip . we have that Om > 0. 



4.2 Connecting the quadratic form to the Sobolev norm 

The benefit of Lemma 14.31 is that we can use it to obtain useful local versions of the form 
{Cu, w)l2(m)- III particular, we consider, for C^n = Q(A), coefficients oq, . . . ,am as guaranteed 
by Lemmag^l When Q C M, the form (EJlo ^)L2(fi) = J2]Lo is the 

local version of {Cu,v)i^^(j^y Indeed, we have 



Jn 

where /3(-,-)x : C~ x ^ M is the bilinear form l3{u,v)x = Y^f=Qaj{V^u,V^v)x. Clearly 
for 

a := max \aj\ and a= max \aj\ (4.1) 

we have 

m—l m 

^{V^u,V^u)x. (4.2) 

j=0 j=0 
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If we integrate over a region Q C M, we obtain 

Now if u vanishes on a sufficiently dense set Xq C 0,, then a corollary of the "zeros" 
estimate Theorem lA.lH given in Section [Aj will imply that Oml^^l^m^j^^ — a||u||^„_i^^^ > 

e||u||^m(Q-) where e depends on am, a, properties of Xq and the boundary of J7, but noth- 
ing else. The two most important types of subset fi, for our purposes, are annuli a and 
complements of balls b^. 

Annuli. In this case we consider an annulus a = B{pQ,R) \ B{pQ,R — t) with outer 
radius R < ^rm and we apply Corollary IA.16I to a function u that vanishes on a set X 
satisfying h = h{X,a.) < jt. If, in addition, h < y/j^, we have < A/i^|n|^mj^^j < 

2^\'^\w^{a.) simultaneously, ^lkll^m(a) < ^[""In^^Ca)' ^° "■"^\'^\w^{n)~^\\'^\\w^-^{n) - 
^\\'"'\\w:^{n) follows and 

2 / / / ttII„,I|2 



— ll^llH-2-(a) < J Ku,u)dn<a\\u\\^^(^^y (4.3) 

(Note that h must be chosen to be less than y^^^, as well as 7^). 

Complements of balls We consider the punctured manifold b'^ = M \ h{pQ,R) with 
outer radius R < ^im- We apply CoroUarv IA.17I to a function u that vanishes on a set X 
satisfying h = h{X, a) < jt. By picking h < y^^^, 

^\Mw^(h(p,rr) < ^^'^'^ - "ll^llH/r(b(P,r)-)- (^"^^ 

follows. (Note that in this case, h must be less than /iq and but that it can be chosen 

independently of r. In Lemma 15.11 we refer to this critical value, the minimum of ho and 
^/^, as Ho.) 



5 The Lagrange function 

We wish to uniformly bound the Lagrange function Xii^) and establish its rate of decay as x 
moves away from its center ^. There are two cases that we will consider. 

The first is the special case that involves interpolation by a polyharmonic kernel km (cf. 
Definition 13. 2p that is conditionally positive definite with respect to a space 11 annihilated 
by the operator Cm- This case is significant because the rate of decay is exponential (cf. 
Theorem 15. Sp . It includes the restricted surface splines on S'^ discussed in Example [21 for d 
even. 

The second case is the general one, where we do not assume any annihilation properties 
concerning the space 11 that is to be reproduced. This case includes the surface splines in odd 
dimensions. The decay rate in this case is algebraic, rather than exponential. 



17 



These results are similar to ones for the case of a lattice in M"^ [5]. The restricted surface 
splines defined in ()3.2p have Lagrange functions that decay exponentially, for d even, but only 
algebraically for d odd. For d odd, the lattice case has an additional family of polyharmonic 
splines with exponential decay. We conjecture that this exponential decay holds for odd- 
dimensional spheres, and that we have obtained only algebraic decay is simply an artifact of 
the proof. 



Notation 


Constant 


Introduce in ... 


(I'm 


(Positive) lead coefficient of the polynomial Q{x) 


Lemma 14.31 


a 


maximum coefficient of Q{x) 
(in absolute value) 




a 


maximum coefficient of Q{x) — UmX™' 




Cq 


£i(M"^) norm of eigenvalues of Cm\i\_j 


(I53D 


I'M 


injectivity radius 


Section [2] 


ri,r2 


constants of metric equivalence from Exp 


(lOD 


Cl,C2 


constants of metric equivalence for Sobolev spaces 


Lemma 12.21 


A 


constant for zeros lemma for annuli 


(1A.20P 


ho 


threshold h level for the zeros lemma 


(A.19) 


Ho 


threshold h level for results of 15.11 


Lemma 15.11 


Hi 


threshold h level for results of 15.21 


Lemma 15.41 



Table 2: Constants frequently used in Section [5j The first four constants are related to the 
elliptic operator Cm = <3(A). The final seven are geometric constants depending on M. 



5.1 Cm annihilates 11^ 

We first consider the special case where km satisfies p.Sp . with an operator Cm = Q{^) for 
which km{j) = iQ{j))^^ > ioT j ^ J' and Cm4>j = for j € J^. In other words, km is 
conditionally positive definite with respect to XIj-, and Cm^j = 0. This the case for Example 
[2] for surface splines on even dimensional spheres. 

In this case, the native space seminorm (jS.ip is precisely the quadratic form derived from 
the operator, namely 

The more general case is considered in the next section, although the basic elements are 
present here. 

We begin by observing that if H is sufficiently dense, with h < min(/io, then by 

(|4.4p it is possible to estimate the norm of the Lagrange function by comparing it to a bump 
(j)^ with (p^ o Expg(x) = (j{\x\/q). We note that this bump is 1 at ^ and vanishes on the rest 
of H, thus it interpolates xs, on ^-nd has a smaller native space seminorm. 

^Wxdw^iM) < ^ i'^?iL,J - ^WkWwriM) < Caq'^-^"'. (5.1) 
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The final inequality follows from Lemma [221 and a direct computation of ||cj(| • |/9)||vi/2m(-]gd). 

The main result, the near-exponential decay of the Lagrange functions, now is a conse- 
quence of an argument developed in [16j but given here in a somewhat different, streamlined 
form. First we prove a lemma showing that a fraction of the seminorm of the Lagrange func- 
tion taken over the punctured manifold b(^, ry resides in a narrow annular region around 
the circle dist(j;,^) = r. 

Lemma 5.1. Suppose M is a d- dimensional compact Riemannian manifold satisfying As- 
sumption \4-l\ Suppose further that m > d/2, km satisfies Definition \3.S\ and that Cm annihi- 
lates the space lij- Then there is a constant K > 0, depending only on m and M so that the 
following holds. If H is sufficiently dense, meaning that 



h < Hq := min ^/iq, 

and if a = h{p, r) \ b(p, r — t) is an annulus of outer radius r < im (^nd sufficient width t, so 
that Ah/ ho < t, then the Lagrange functions for interpolation by km satisfy 

\\Xi\\w:^ih{p,r-t)'^) ^ ^\\Xi\\w^ibip,r)\b{p,r-t))- 

Proof. Since xs, minimizes the native space seminorm we have [xc]fc,„ M — I'^^CX?]!™ M 
function 0g equaling 1 at If 0^ is a cut-off, equaling 1 in the ball b = h(p, r — t) and 
vanishing outside of the ball b U o, then 

By (gaD and (gS]) 




The result follows with K = 4:aK'/am, where the constant K' is introduced in Lemma [521 
which we prove below. □ 

Lemma 5.2. Assume the manifold M, the kernel km, the set of centers S and the annulus 
C M satisfy the conditions of Lemma \5.1\ If (j)^ is a smooth "bump" function, satisfying 

fl 2t — r\ /\x\ 2t — r\ 
4>^ o Exp5(x) = a f -dist(Exp^(x),Exp^(0)) —\ = ^ ( ^ + ^^'^^ 

with a : IR+ — t- ]R-|_ a C°°, non-increasing cutoff function equaling 1 on [0, l\and on [2,oo), 
then 

Whxdw^ia) <K'\\x^\\w--{a) 
where K' depends only on M, m and the choice of cutoff a. 
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Proof. We follow the proof of [16\ Lemma 4.3]. Let X(.{^) = X(, ° Exp^. By using the metric 



equivalence guaranteed by Lemma 12.21 we can estimate 



|2 < „2 

\ Co 

'2 ~ 



E 

|a|<m 



x\ 2t-r\^, 



2 

dx 



by 



< c: 



< 



j=0 

)2 m 
j=0 



B{0,r)\B{0,r-t) 



m\wi(a) 



\»\=j 



and K' = C(^)^A ^J1q(/io/4)^(™~-'\ The second inequality follows from the product rule, 
and C is a constant depending only on m, d and a. The third inequality is Lemma 12.21 again . 
and the final inequality is the zeros lemma for annuli, Corollary IA.16[ □ 

At this point, we can follow the example of [161, Section 4] 

Theorem 5.3. Suppose that M is a compact d- dimensional Riemannian manifold satisfying 
Assumption \4- 1\ Suppose further that m > d/2 and that km satisfies Definition \3.2\ and that 
Cm annihilates the space Hj. There exist positive constants ho, v and C, depending only on 
m, M and the operator Cm so that if the set of centers H is quasiuniform with mesh ratio p 
and has density h < Hq then the Lagrange functions for interpolation by km satisfy 



\xdx)\ < Cp'—^/^exp (-^min((i(x,0,rM)) • 



(5.3) 



Furthermore, for any < e < 1, there is a constant C depending only on m,Wl, p and e, so 
that the Lagrange functions satisfy 



\xd^)-xdy)\<c 



d{x,y) 



(5.4) 



Proof. Set t = Ah/ho =: 7/1, and note that for t < r < im, Lemma |5 . 1 1 implies that 

IIXcllH/2'"(b(5,r)'=) - ^IIX5llw'2'"(b(^,r-t)=)' 

with e = {K — V)/K. Letting n := [r/tj, we have 



< f:-lp(l°§')'"/*llvJ|2 

< Ce'''''/^Q^'^-'^ 



^ -1 -urlh\\ ||2 



with V := — 7loge. Since e = < 1, it follows that > 0. The final inequality follows 
from (EH). 
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The bound ()5.3p follows from the observation that X({x) can be estimated using The- 
orem [AyTIl The intersection h{x,t) n b(^, dist(x, ^))'^ is contained the ball h{x,R), since 
t < R < irM- Because geodesic spheres are smooth hypersurfaces whose intersection is 
nontangential, The intersection h{x,t) D b(^, dist(j;, i^))'^ is a Lipschitz domain contained in 
h{x,R). Moreover, h < ^t. Thus, Theorem lA. Ill applies, giving us 

for h < 'jt. Similarly, the estimate in (j5.4p follows from Corollarv lA.151 □ 



5.2 General Case 

In this case, the native space seminorm (j3.ip and the quadratic form induced by the operator, 
differ by some low order terms: 

Because of the orthonormality of 99^, we have \{u,(pj)\^ < W^W^^^^^^y Setting 

Cq:=Y,\Q{W\ (5-5) 

(this is the ^i(]R^) norm of the spectrum of the operator Cm restricted to Hj-) we note that, 
by Corollary IA.13t if u vanishes on a sufficiently dense set, then the lower order terms are 
controlled 

Indeed it follows that ^||w||pt/m < 111^^11^^ j when h is chosen small enough that Cqh?^ < 

This allows us to provide a basic estimate for the Lagrange function, similar to ()5.ip . In 
this case 



a 



-fWxdk-m ^ < miLj < ia + CQh'n\m\w,-iM) < Caq'-'-^- (5.6) 

Lemma 5.4. Suppose M is a d- dimensional compact Riemmanian manifold satisfying As- 
sumption \4-l\ Suppose further that m > d/2 and that km satisfies Definition Then 
there is a constant K > 0, depending only on m and M so that the following holds. If H is 
sufficiently dense, meaning that 



h < Hi := min ( Hq, 2™/ ) = min ( Hq, 



2Aa' ySCqJ 

and if a = h{p, r) \ h(p, r — t) is an annulus of outer radius r < tm, satisfying, in addition. 
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(for a constant Cq depending only on m, M, k„i, and J which we define in the proof) and 
sufficient width t, so that h < 'jt, then the Lagrange functions for interpolation by km with 
auxiliary space Hj satisfy 

IIXcllw'2'"(b(p,r-t)=) ^ ^\\Xi\\w:p{b{p,r)\b{p,r-t))- 

Proof. Since xs, minimizes the native space seminorm, we have 

2 



for a cut-off (j)^ equahng 1 in the bah b = h(p, r — t) and vanishing outside of the bah 
b(p, r) = bU a. Using the sum of squares factorization \A\'^ — \B\'^ = ^[{A — B){A + B)~\ , we 
may write 



(f)^{x)x^{x)(pj{x)dn{x) 



X^{x)(pj{x)dfi{x] 



^) - l)x5(a;)v'j(a;)d/i(a;) x / {(t>^{x) + l)x^{x)(pj{x)dn{x) 



{M^))xd^)fj{x)dfi{x)- / xd^)'Pj{x)dfi{x 



X ( / iM^) + '^)xd^)^j{x)dii{x) 



The second factor can be bounded by using Coroharv lA.lSl along with the cutoff function 
being bounded by 1 and \\(pj\\2 = 1: 



+ l)xg(x)(/?j(x)d/i(x) 



< 2||xdU,(M) < 2A/i'"||xdll4',^ 



•(M) 



To bound the first factor, start by using Corollary IA.16I and Lemma 15.21 to obtain 

1/2 



{4>d^))xd^)'Pj{x)dn{x] 



< Ah"'U^Xi\\wria) < ^K'h^\\Xi\\wrr.^,y 



Next, from Corollarv IA.17I we have that 



X^{x)fj{x)dfi{x] 



1/2 



< ( I^Jxil'df^ix)] <Ah"-\\x^\\wi^i,cy 



So the first factor is bounded by 

AK'h^\\xdw,-ia)+^h^' 



vi/-m((,c) < A{K' + l)h'^\\x(\\w.^(b<=), 
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and the product itself is bounded by 
Putting these bounds together gives us 



, where C = 2k^{K' + 1). 



(j)^{x)x^{x)^pj{x)(l^l{x 
< C' CQh^'^\\xii\\w^{M)\\X£\\w^{b'') 



X^{x)ipj{x)dii{x] 



Thus for h sufficiently small, say for C CQh?"''\\x£^\\w:!p{M) ^ ^llx^llvKj^Cb'^)) which follows 



by taking 
we have 



• = 



(/)5(x)xc(x)(/?j(x)d//(x) 



X£_{x)ipj{x)dn{x) 



< 



< 



""in II II 2 {K 7\ 

(5.8) 



We note from (|4.4p that 

-^\\xdw:^ib'') 

and by subtracting right and left sides of (j5.7p from the left and right sides of (jS.Sp the lemma 
follows since then 



,^{b-) 

2 



X^{x)(pj{x)dfi{x] 



<P(ix)xd^)Vj{x)dfi{x) 

< [kx^]lm,b- = [kx(]l^,a < ^K'Wxdw^ia) 

where the last inequality follows from (j4.3p . The result follows with K = SaK' /am- □ 

We are now ready for the full result. 

Theorem 5.5. Suppose that M is a compact d- dimensional Riemannian manifold satisfying 
Assumption \4-l\ Suppose further that m > d/2 and that km satisfies Definition \3.S\ There 
exist positive constants Hq, v and C , depending only on m, M and the operator Cm so that 
if the set of centers H is quasiuniform with mesh ratio p and has density h < Hi then the 
Lagrange functions for interpolation by km with auxiliary space H j satisfy 



\xdx)\ < Cp"-'^/2max(exp(-^min(d(x,C),rM)) ,h^'^ 



(5.9) 



Furthermore, for any < e < 1 for which m > d/2, there is a constant C depending only on 
m, M, p and e, so that the Lagrange functions satisfy 

'd{x,yy 



\xi{x)-xdy)\<c 



(5.10) 
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Proof. Let rg be the smallest radius r so that ||x5llw2'"(b(p,r)':) < C'o/i^"*||x^|| W2™(M)- Since 
^ \\xs.\\w:^{h{p,ry) is decreasing, tq < diam(M). Assume without loss that tq < rM, since 
otherwise the proof proceeds exactly as in Theorem 15.31 

Set t = h/'j, and note that for t < r < tq, Lemma 15.41 implies that 

IIX5llH/2'"(b(5,»-)'=) - ^\\X^\\w^{h{i,r-t)'') 

with e = {K-l)/K. 

As in the proof of Theorem 15.31 

II ||2 ^ —l—vr/h\\ ||2 ^ /-I —VT Ih 2m— d 

\\Xi\\w^ib{^,rY) ^ ^ ^ 1 1 1 1 1^2™ (M) 1 

Where we have set v := — 7loge. Since e = < 1, it follows that u > 0. The last 

inequality follows from (j5.6p . On the other hand, for r > r^, we have that ||xcllH/'"{b{p,r)=) ^ 
Coh^'^WxdwriM) < CCo'7^"-'^/i^™, by Therefore, 

\\xdwr(Mp,rr) < C<72™-«'max(/i2-,e— Z'^) 

Again, estimate (|5.9p follows from the observation that x^i^) can be estimated by way of 
the zeros lemma: 

\xd^)\ < C'/i'"~'^^^IIX€llw/2"(b(x,t)nb(€,dist{x,0)^) ^ '^^™"'^^^ll^«llH/™(b(e,d(x,o)')' 
for h < 7t. 

Similarly, estimate (jS.lOp follows from Corollary lA. 151 □ 



5.3 Implications for Interpolation and Approximation 

At this point, we are able to state three important corollaries to Theorem 15.51 that satisfac- 
torily answer the questions concerning bases and approximation properties of Vx discussed 
in Section [TJ These results were previously obtained in \16\ 115] for a class of Sobolev kernels. 
Here, we get them for a much broader, computationally implementable class of kernels. Our 
first result is that the Lebesgue constant for interpolation is uniformly bounded. 

Theorem 5.6 (Lebesgue Constant). LetWl be a compact Riemannian manifold of dimension 
d satisfying Assumption \4-l\ Suppose further that m > d/2 and that km satisfies Definition 
\3.2[ For a quasi-uniform set H C M, with mesh ratio h/q < p, if h < Hq, then the Lebesgue 
constant, L = sup^g^ ^^g^ |xg(a)|, associated with km and J is bounded by a constant 
depending only on m, p and M. 

Proof. Fix X. Using Theorem 15. 5( we estimate the sum as 
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The first sum can be treated exactly as in |16^ Theorem 4.6], and is bounded independently 
of h. The second sum, II, can be estimated using the fact that < Cq~'^, with a constant 
C = fi{M)/a{M), where a(M) := inf^gM info<r<rM '^~V(t»(x, r)). Thus 

which is bounded since 2m > d (indeed it vanishes as /i — )• 0). □ 
The next consequence is the Lp stability of the Lagrange basis. To this end, we define 

S{km, J, ^) = \ Y1 6 + P I P e and ^ A^q{0 = for all qGUj 

and use this notation in lieu of Vx used in the introduction. 

Theorem 5.7 (Stability of Lagrange Basis). Under the assumptions of Theorem 15.51 there 
exist constants < ci < C2, depending only on km, J , M and p so that 

ClPp,-ll^(H) < I|s||Lp(M) < C2pp,.||^p(H)- 

holds for all s = ^^X^ ^ S{kfn, , with normalized coefficients Ap^^ := q'^/^A^. 

Proof. When Cm annihilates lij, this is a direct consequence of the pointwise estimates 
obtained in Theorem 15. 3[ We observe that the following three conditions hold: 

1. The basis (x^)^^^ is a Lagrange basis. 

2. The basis has decay |xc(2;)l ^ C/o™'"'^/^ exp (— ^ min((i(x, ^), rMi)) . 

3. The basis has the equicontinuity condition 1x^(2;) — X£,{y)\ ^ C2 



dist(a::,y) 
9 



Thus the result [151 Theorem 3.10] applies. 

In the general case, the result still holds, despite the fact that item 2 may fail. Le., the 
Lagrange functions may decay more slowly than the basis functions considered in [TS], and a 
minor modification is required to apply of [151 Theorem 3.10]. 

The upper bound \\s\\]^ < C2||Ap^.||£ follows directly from the estimate (j5.9p . Indeed, 
the case p = cxd, is none other than the Lebesgue constant estimate Theorem 15.61 while the 
p = 1 case follows by the uniform bound on 

llxdii < C'P""'^/^ (Ch'^ + vol(M) + e-^'^/^^y^ < Cp^^'^l'^q^. 

The case 1 < y> < 00 follows by interpolation. 

To handle the lower bound, we utilize functions c/)^, defined in a similar way as in (|5.2p . 
satisfying o Exp^ = o", with 



a{x) 



1 \x\ < ro 



h-^'^e h ro < \x\ < im or a{x) = , , (5.11) 




\x\ < ro 
ro < \x\ 
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and with threshold value tq := —^hlogh (the second definition is chosen if rM < tq). It 
follows that 

= Xii4>i + X€(i - (P^) =■■ 55 + 

and satisfies items 1-3 above (in particular, item 3 follows since (f)^ is bounded and Lip(l), 
with Lipschitz constant I'/h) and |15l Theorem 3.10] applies. In particular, there is ci > so 
that II E^eH^e^cllp > ciPp,.||£p(H)- 

On the other hand, \b^ix)\ < /i^"", implies that || Z^^gs^^^ellp ^ Cp'^lWAlip^h'^"''''' since 



/ |^A^6g(x)|dx < p||^j(H)^(M)/i2"^ and max\"^A^b^{x)\<Cp'^\\A\ 
Jm 

Thus, for s = X^geH 



1,2m— d 



1/p 

'Cl 



where we have used the inequality | ^^g^ ^^5? 1^ — ^^^^ (l X^geH ^^X^l^ + I ~ S^eH 

□ 

Our final consequence treats the Lp stability of the L2 projector. This was a primary 
goal of [15j, and, in light of Theorem 15.51 we can produce a similar result here, with a minor 
modification to handle the slower decay of the Lagrange functions. 

Let V : — )• S{km, ,'^) be a basis "synthesis operator" V : (^^)gg= 1— )■ X^^gs^^^^^, 
for a basis (u^)^g= of S{km-,J Likewise, let V* : Li(M) — t- be its formal adjoint 
V* : f ^ {{f,v/:))^^^. The L2 projector is then Th = V{V*V)-^V* : Li(M) ^ S{km,J,E), 
in the sense that when / G L2(M), Ts/ is the best L2 approximant to / from S{km, J ,'^)- 

The L2 norm of this projector is 1 (it being an orthogonal projector), while the Lp and 
L'p norms are equal, because it is self-adjoint. Thus, to estimate its Lp operator norm (1 < 
p < 00), it suffices to estimate its L^o norm. 

Theorem 5.8. Under the assumptions of Theorem \5.5\ for all 1 < p < 00, the Lp operator 
norm of the L2 projector is bounded by a constant depending only on M, p k^ and J . 

Proof. When Cm annihilates 11 j , Theorem 15.31 and Theorem 15.71 satisfy the conditions of |151 
Theorem 5.1] (the Lagrange basis is stable and rapidly decaying), and the result follows. 

In the general case, we cannot directly apply this theorem, because the basis does not decay 
rapidly enough. We take as our basis = X£,,2 ■= Q~'^^'^X^7 the L2 normalized Lagrange basis. 
It follows from Theorem [521 that \\V\\e^(^E)^L^{M) < C2q~'^/'^ and ||T^*||loo(m)^£oo(h) < C2q'^^'^ ■ 
Thus, to estimate the L^o operator norm of (and thereby all other Lp norms), it suffices 
to estimate the ^00 (^) ^00 (^) norm of the inverse Gram matrix {y*V)~^ . 

We make the split g^ = Xg^g and b^ = X£_0- — (^e) with 0^ o Exp^ = a defined as in (j5.1ip . 
And note that X5,2 = Xc,2(/'c + X?,2(l - (t>^) =■ 5^,2 + ^5,2, 



26 



It follows that V*V = G + B, with G^^,^ = (5^,2; fl'c.2)-L2(M)- 

The functions {g^) are a Lagrange basis, in the sense that 55(C) = (although they span 
a different space than S{km, J, and, as observed in the proof of Theorem 15.71 they are Lp 
stable. They also satisfy the decay conditions of [15, Proposition 4.1], and by applying this 
result we see that ||G~^||oo is bounded by a constant. 

On the other hand, IS^,^ < | (5^,2, &c,2)l + K&5,2, 9C,2)| + K&5,2, &C,2)I < C/i^*^, and \\B\\oo < 
max^Ec ^ Cp'^h'^"'-'^. The theorem follows by noting that {V*V) = G(Id + G'^B), 
and, hence, \\{V*Vy^\\oo < \\G-^\\oo (1 + o{h)) . 

□ 

5.4 Spheres and S0{3) 

We now explore some further consequences of the results of the previous section. We will 
shortly see that Theorems 15.61 and 15.81 imply that 1= and Ts are near-best. In some important 
cases, we can then use these projectors to observe precise rates of convergence for interpolation 
and least squares minimization, better rates than were previously known. 

For the kernels considered in section 13.21 theoretical approximation results are known 
in some special cases, including spheres and SO (3). The difficulty is that these results are 
often not practical, because they are derived from approximation schemes that are difficult 
to implement. The good news is that the stability of the schemes /= and Ts imply that these 
operators, which are associated with practical schemes, inherit the same convergence rates. 
Indeed, for a normed linear space Y and a bounded projector P :Y ^ Y, one has for f £Y, 

\\f-Pf\\= inf IIZ-s + Ps-P/ll < (l+||P||)dist(/,ranP). (5.12) 

s£ran_P 

This fundamental observation is known as a Lebesgue inequality, and we employ it with 
P = 1= and Y = C(M) as well as with P = T'^ and Y = Lp(M). In recent years, a concerted 
effort has been undertakerB to understand the general Lp convergence rates (i.e., the behavior 
of dist{f , S {km, J ,'^))p as H becomes dense in M) of certain well-known kernels in terms 
of smoothness assumptions on the target function / and on the density of the point-set S, 
measured by the fill distance h. 

To measure smoothness of the target function, we make use of the classical (Sobolev, 
Besov) smoothness spaces introduced in Definition 12.11 and Definition 12.31 with the exception 
that for approximation in L^o, we make the (usual) replacement of C"^"^ for (but using 

the same norm). As a shorthand, we capture the smoothness spaces we use by means of a 
common notation, Wp(M). For m > d/2 denote the space Wp(M) by 

• >V^(M) = C2'"(M), when p = 00 and s = 2m 

• W*(M) = W^'^iM) when 1 < p < 00 and s = 2m 

*This stands in contrast to the more classical, mainstream theory of kernel approximation, where approx- 
imation properties are investigated and understood only for functions coming from the reproducing kernel 
(semi-)Hilbert space associated with a conditionally positive definite kernel. Obtaining an understanding out- 
side of this context has generally required indirect, theoretical approximation schemes, and it has not been 
obvious, until now, that such results would have practical consequences. 
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• yV^(M) = B^^^{M) when 1 < p < oo and < s < 2m. 

Corollary 5.9. For M = S'^, and m > d/2 the surface splines introduced in Example 3: 
km{x,y) = (f)s{x-,y) satisfy the following. There is a constant C so that, for a sufficiently 
dense set E C S'^, and for f e . 

1. For f G W^, 11/5/ - /II < Ch'WfWw^ 

2. Forl<p<^ and for few;, HTh/ - /||p < /i^||/||w|„. 

Proof. As with the Sobolev kernels, (j)s is of the form Gp + -0 * Gp as considered in [22], the 
result follows direct from |22[ Theorem 6.8]. Alternatively, it follows from \\4\ Theorem 6.1], 
which treats kernels on the sphere of the type in Definition l3.2[ The Besov space result follows 
from [2?, Corollary 6.13] or [HI Corollary 6.2]. □ 

Corollary 5.10. ForWi = SO{3), andm > 2 The surface splines, k^, introduced in Example 
4 satisfy the following. There is a constant C so that, for a sufficiently dense set H C S0{2>), 
and for f G Wp . 

1. For f G W^, ||/h/ - /II < C/ill/llw^ 

£ Forl<p<oo and for few;, ^f - f\\p < h'\\f\\wg^. 

Proof. This follows from Theorem 9] for the case of full smoothness and from [17', Theorem 
12] when < s < 2m. □ 



A A Zeros Lemma for Lipschitz Domains on Manifolds 

Results concerning Sobolev bounds on functions with many zeros are known for Lipschitz 
domains in \2,h\ [26] . Our aim is to extend these results to certain Lipschitz domains on 
manifolds. Before we do that, however, we will need to improve the results in [25|[26]. 

A.l Lipschitz domains in 

Consider a domain 17 C M"^ that is bounded, has a Lipschitz boundary, and satisfies an interior 
cone condition, where the cone Cq has a maximum radius and apertur^ 299. Of course, 
the cone condition will be obeyed if we use any radius < R < Rq. The theorem that we will 
give below requires covering with certain star-shaped domains. 

We will say that a domain T> is star shaped with respect to a ball B[xc,r) := {x G 
W^: \x — Xc\ < r} if, for every x G P, the closed convex hull of {x} U i? is contained in P [H 
Chapter 4]. For T> bounded, there is a measure of how close to spherical P is; namely, the 
chunkiness parameter 7 [H Definition 4.2.16]. This is defined as the ratio of diD to the radius 
of the largest ball relative to which T> is star shaped. When I? is a sphere, 7 = 2. If there is a 
ball B(xc,R) 5 P, then r < d-p < 2R and 7 < Finally, such domains satisfy an interior 
cone condition and certain Sobolev bounds, which are stated in the next two propositions. 

^Aperture here is the angle across the cone, 2(p in this case. In optics, aperture would be (p. 
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Proposition A.l ( \25\ Proposition 2.1]). IfD is bounded, star shaped with respect to B{xc,r) 
and contained in B(xc,R), then every x ^T) is the vertex of a cone Cx> C P having radius r 
and aperture 6 := 2 arcsin (5^) ■ 

Tliis proposition also implies that the chunkiness parameter for D is bounded in terms of 
the aperture: 

7 < — = csc(6'/2). 
r 

Proposition A. 2 {fW, Proposition 3.5]). Let V C be as above, m G N and p G M, 
1 < p < 00. Assume m > d/p when p > 1, and m > d, for p = 1. If u ^ W^{T>) satisfies 
u\x = 0, where X = {xi, . . . ,xn} C T> and if h = hx < je^p:^; then 

IhllL^(D) < C'm,d,p7'^c?o~'^'^Vlw'p-(©)- (A.2) 

Our next task is to obtain Sobolev bounds for the domain 17 C M'^ that are similar those 
in (jA.ip . The idea is to cover Q with star-shaped domains. To do that, we will use a 
construction due to Duchon [8J . With Rq , 293 being the radius and aperture for the cone Cq , 
and < < Rq, let 

sinfw) , ^ f 2r „ . 

„„J rr J + ^ '77a. 



r := 2RF((p), where F(m) := — and Tr:=<te -^ll" : B(t,r) CU} . (A.3) 

4(l + sm(vp))' \ Vd J 

For t G Tr, let T>t be the set of all x G such that the closed convex hull of {x} U B[t,r) 
is contained in Q, n B(t,R). From |25l Lemma 2.11], we have that each Dt is star shaped 
with respect to the ball B{t, r), and satisfies B(t, r) CDt C. 0,0 B{t, R), dx>t < 2R. Because 
2R/r = 1/F{lp), the aperture for Cx>t is 

e = 2arcsin(l/F(v3)), 

and the chunkiness parameter 7^ for T>t is uniformly bounded: 

2 < 7* < — = (A.4) 

We also have that Q. = UteT, '^t^ t^at < Cd^o\{Vt){F{if)R)-'^ , and that vol(A) < CdR'^. 

The integer- valued simple function "Ylit^Tr XB{t,R){^) is the number of B{t, i?)'s that contain 
X. This is easily bounded above by M^^ip, maximum number of such balls intersecting a fixed 
one, say B{0,R). A little geometry shows that 

M{d, if) < {2R/r + l)'^ < 2'^/{F{ipf 

Note that the existence of Ma^i^ implies that for any function / in Li($7) we have 

I \f{x)\dx= I Y.^vAx)\f{x)\dx<Md,^ [ \f{x)\dx<i2''/{Fiip)^) [ \f{x)\dx. 
^ JVt Jo, I Jn Jn 

(A.5) 
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Lemma A. 3. Suppose that h = hx,n satisfies h < ^^F{(p)'^, then I! A . 1\) and liA.l^) hold 
uniformly in t for Df, provided 'jt and dx>t are replaced by 1/F{ip) and 2R, respectively. 

Proof. The mesh norm for Vl satisfies 

h < m-'^{2RF{ip)) X {F{^)/IQ) < r, 



since F{(p) < 1. It foUows that B{xc, r) n X ^ 0, and so 2?j n X contains at least one point of 
X. From this we have that hx>tnx ^ h. The lemma then follows from the bound on h being 
less than the one required in Proposition IA.2[ □ 



We wish to prove the following result, which differs from an earlier result in \25\ Theo- 
rem 2.12] in that it applies to cases in which the index A; < m — 1, as opposed to k < m — n/p. 

Theorem A. 4 (Euclidean Case). Suppose that Q is a Lipschitz domain obeying a cone 
condition, where the cone Cn has radius Rq and aperture 2(p. Let k, m, and p be as in 
Proposition \A.S\ and and let X be a discrete set with mesh norm h satisfying 

H < ^,Fi^f. (A.6) 



If u ^ W^{^) satisfies u\x = 0, then 

p(^^ym+d+d/p 



and 

Proof. Given h, choose R = Sm?h/ F{ip)'^ < Rq. Applying Lemma I A. 31 and Proposition IA.2I 
to the domain T>t then results in the bound 



FiifY 



We will follow the proof in [25, Theorem 2.12]. Summing over t on both sides of the previous 
inequality, using $7 = V}{Dt and applying ()A.5p . we have that 



from which ()A.7p is immediate. The bound on 11^11^^(0) follows similarly. □ 
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A. 2 Lipschitz domains in M 

A domain on a smooth, compact Riemannian manifold M satisfies an interior cone condition 
if there is a cone C C M" with center 0, aperture 2(p, and radius R such that, with some 
orientation of C, Exp^ : C — )• Cp C 0. That is, the image of the fixed cone C is a geodesic 
cone Cp contained in Q. In addition, Q satisfies the uniform cone condition if, for every 
Po G dCl and some orientation of C, Expp(C \ {0}) C $7 for ah p £ h{pQ,r) n $7. Finally, 17 is 
said to be locally strongly Lipschitz [231 119j if for every pQ G there is a local chart (U,^), 
i) -.U M",with V'(po) = 0, a Lipschitz function A : M"-^ M, with A(0) = 0, and an e > 
such that 

i^{U r\Q) = {(x, A(x') + t) : < t < e, x' G M"-\ < e}. 

Our approach to a manifold analogue of Theorem IA.4I is to employ a set of points for M 
that are similar to those described in (IA.3h . The set that we need is described and studied 
in [131 §3]- Let e > 0. There exists an ordered set of points {pi, . . . -.pn} C M such that 
the uj*^]^b(pj, e) = M and such that the balls b(pj,e/2) are disjoint. Such a set is called a 
minimal e-net in Mllfl It has the following two important properties: First, there is a number 
Ni = iVi(e,M) for which N < Ni. Second, there exists an integer N2 = iV2(M) > 1 such that 
for any p G M the ball h{p, e) intersects at most N2 of the balls h{pj,e). It is remarkable that 
N2 is independent of e and, in fact, depends only on general properties of M itself. We will 
need a slightly stronger version of this result. 

Lemma A. 5. Let {pi, . . . ,pn} be a minimal e-set, p G M, and let 1 < q. Suppose e < dyi/a, 
where dm is the diameter o/M. Then the cardinality s := #{pj- h{p,ae) n h{pj,s) 7^ 0} < 

3{d-l) , 

(4a + l)'^eW". 

Proof. The argument used in |13^ Lemma 3.3] gives, mutatis mutandis, 

ff"^'^^' smh'^-\J\^t)dt , , rr^ , 

J^'smh'^'^{y^\t)dt 

where (d— 1)k is a lower bound on the Ricci curvature of M. Making use of 1 < sinh(a;)/x < e^, 
we see that 

d-^x'^< r^^-^dtK r smh^-\t)dt<d-^x'^e^'^-^^'', 
Jo Jo 

and consequently that 

3(d-l) J 

Hia, e, k) < (4a + i)'^e('^-i)(2a+i)./^ < ^ ^^d^ ^ '^m^ 

which completes the proof. □ 

''An e-net is a set of points X = {pi,...,pjv} for which |Jb(pj,£) covers M - in other words, for which 
h{X,M.) < e. Likewise, a minimal e-net is quasiuniform, with separation distance q > e/2 and mesh ratio 
h/q < 2. 
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Lemma A. 6. Let i? < tm/S, G (0, 7r/2] and e = 2^i+sTn(y)"*) • If {pit ■ ■ ,Pn} is an e-set and 
if Cp is a geodesic cone with center p, radius R, and angle then for some j we have that 
h{pj,e) C Cp. 

Proof. We will work in normal coordinates on TpM, where the cone C has vertex at the origin 
and en, = (0, . . . , 1) is chosen to be along the axis of C. The largest Euclidean ball in C has 
radius p = Rsm{ip)/{l + sin{(p)) and center Xc = {R — p)en- It follows that any ball having its 
center a Euclidean distance p/2 from Xc and having its radius less than p/2 is also contained 
in C. Let Pc = Exp(xc). Since the balls h{pj,e), j = 1, . . . , N , cover M, we can find pj such 
that Pc G h{pj,e). 

Let Xj = Expp^(pj). Equation (|2.3p implies that |xc — Xj| < dist{pc,Pj)/ri < e/Li = p/2. 
Now consider the ball h{pj,rip/2). Let q G h{pj,rip/2) and let x = Exp~^(g). Applying 
equation (j2.3p then yields that jx — Xj| < dist(j), q)/Ti < p/2, and consequently that h{pj,r) C 
Cp, with r<p/2 = □ 

Our goal is now to cover with domains analogous to those used in the previous section. 
To that end, let R < Rq, fix 

Fii^sinfoj) ^ „ , . X 

r := — i = 2F(ip)R (A.9) 

2(l + sin(v9)) ^ ^ 

and find a minimal e-net (with e = r) {pi, . . . ,pn} and set Tr := {pj : h(pj,r) C 11}. Because 
obeys a uniform cone condition, with radius Rn and angle (p. Lemma I A. 6 1 implies that 
is nonempty. 

Next, for each pj G T,., let T>j be the set of all p G H h{pj,R) such that the geodesic 
convex hull of {p} U b(pj,r) - i.e., the set comprising all points on every geodesic connecting 
p to a point in b(pj, r) - is contained in n h{pj, R). Again by Lemma \A.6\ for every p G Q 
there is a pj G Tr such that the geodesic cone Cp contains h{pj,r). Since this cone also 
contains the geodesic convex hull of {p} U h{pj,r), it follows that p G Vj and, hence, that 

We claim that the domain Dj := Exp~^(Pj) is star shaped with respect to the Euclidean 
ball b(Exp~,"^(pj), r/r2). To show this, we will need the following lemma. 

Lemma A. 7. Let p £ M, u, v G TpM satisfy \u\p = |v|p = 1, a := arccos((u, v)) G (0,7r]. 
If Pp = Expp(pu), so that p = dist{p, Pp) < rM/3, then the geodesic distance r from pp to the 
ray along v satisfies 

ripsin(min(a, 7r/2)) < r < r2/0sin(min(Q, 7r/2)) . 

Proof. Consider the sector in span(u, v) formed by tu + sv, where s,t > 0. We will work 
in normal coordinates based at p. The minimum geodesic distance r from pp to geodesic 
Expp(sv) occurs at a point Expp(tv). In addition, the minimum Euclidean distance r' from 
pu = Exp~^ [pp) to the ray will occur at another point, t'v, where v is perpendicular to t'v— pu, 
in the Euclidean sense. These facts imply that r = dist(pp,tv) < dist(pp, Expp(t'v)) < 
r2|pu — t'v|cuci- Using a little trigonometry, together with t'v — pu and v being perpendicular, 
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we see that |pu — t'v|euci = psin(a) when a < 7r/2, and that |pu — t'v|euci = |/Ou|cuci = P when 
a > Tr/2. In a similar way, we have r = dist(p, Expp(tv)) > Fil/ou — tv|cuci > Tilpu — t'v|euci = 
ri/9sin(min(a, 7r/2)) . Combining the inequalities completes the proof. □ 

There is corollary to the lemma that will be useful for smooth surfaces, in particular balls and 
annuli. We state and prove it now, although it will only become useful after the zeros result 
Theorem lA.lli 

Corollary A. 8. Let q £ dil. and suppose there is a ball h{p,p), p < rM/3, such that h{p,p) C 
n and that d{q,p) = p. Then, the geodesic cone Cq, with vertex q, axis along the geodesic 
joining q to p, radius Tip/2 and angle = arcsin(2f^) satisfies Cq \ {q} C il. 

Proof. Let q' be a point on the lateral side of the cone that is p away from p - denote it in 
coordinates around qhj q' = sv. With this, we identify two triangles. 

o P 




Figure 1: In this figure, p is the center of the ball of radius /), g is a point on the boundary, 
q' is a point simultaneously on the boundary of the ball and on the side of the cone and p' is 
the nearest point on the ray Expg(tv) to q' . 

The first triangle has corners p = Exp^(/)u), q' = Exp^(sv) and a point on the ray Exp^(tu) 
with < t < /9 ( Lemma A. 7 guarantees that the vertex at p is acute, since s < Tip/2). Let 
us denote the third corner of this triangle by p' = Expg(t'u) Note that p' is a distance of t' 
from q and a distance of p — i' from p. The triangle inequality gives us that 

p < dist(p', q) + {p — t') — > t' < dist(p', q). 

The second triangle we consider has corners g, q' and p' , and Lemma A. 7 ensures that 
dist(g',p') < r2SsinQ! = s/2. So the triangle inequality here gives us that s < dist(p', q')+t' < 
s/2 + t' - so s/2 < t'. 

Combining estimates from both triangles, we see that s/2 <t'< dist(y, q') < s/2, so both 
t' and dist {p',q') are s/2. 
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In other words, the curve from p to q' has the same length as the curve from p to p' to q' . 
By [71 Corollary 3.9, p. 73], any piecewise differentiable curve joining two points - p to p' to 
q' our case - with length less than or equal to any other such curve is a geodesic. Because this 
occurs inside b(p, rfvi), where geodesies do not cross, there can only be one geodesic joining p 
and q' . Since p to p' has to be on the geodesic joining p' to g', and since the length is p, 
and q coincide. □ 

Proposition A. 9. The domain Dj := Exp~,^(Dj) is star shaped with respect to the Euclidean 
ball B(Exp„^ {pj) ,Tir /T2) ■ Also, the chunkiness parameter and diameter for Dj satisfy 



2TlR 



ID, < -p 



£2 

ri 



4(1 + sin(v9)) 
Ti sin((^) 



T.2 

2 



„ , , and dn. < 2R. 
TlF{ip) 



(A.IO) 



Proof. We begin by fixing a point p £ Dj. The geodesic convex hull of {p} U h{pj,r) contains 
a largest cone with vertex p and central axis the geodesic ray connecting p to pj. On this 
cone, whose (lateral) surface consists of geodesies emanating from p, there exists a geodesic 
lying tangent to the sphere dh{pj,r). In other words, we take the cone of largest aperture 2a 
for which all geodesies pass through h{pj,r). 



p 





( ^ 























Figure 2: The largest cone with vertex p and with central axis the geodesic that connects p 
to Pj for which no geodesic lies outside of the ball b(pj,r). The radius r' of the ball centered 
at Pp lying tangent to the cone is greater than 

Prom this two things follow. First, by Lemma IA.71 the distance from p to pj, pj = 
dist(p,pj), the angle a, and the radius r are related by r < r2/5j sin(min(a, 7r/2)) . Second, 
for p < Pj and for a point pp lying a distance of p from p along the central axis, the distance 
r' of Pp to the surface of the cone satisfies: 

r' > ri/jsin(min(a, 7r/2)) > ^ ^ . 

\-2P0 
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It follows the cone contains the ball h{pp, {T2Pj)~^Tirp), which obviously is also contained in 
the convex hull of {p} U h{pj,r). 

Shifting to normal coordinates centered at pj , rather than at p, we see that the geodesic ball 
h{pp, {T2Pj)^^Tirp) contains the (image of the) Euclidean ball B(Exp~^{pp), {T2Pj)^^Tirp). 
Indeed, if \y — Exp~^{pp)\ < (rl/CjO^^Firp then 



dist(Expp,(y),pp) < Tsly - Expp^,^(pp)| < 



A straightforward argument using Euclidean geometry implies that the Euclidean convex hull 
of {Exp~.^(p)} U i?(0, rir/r2) is contained in Dj. Hence, in the Euclidean metric, Dj is star 
shaped with respect to the ball i?(Exp~,-^(pj), rir/r2). Moreover, since Vj C h{pj,R), we 
have that Dj C B(Exp~^{pj), R). Finally, from these facts it is easy to see that the bounds 
in (TOO]) hold. 

□ 

Applying this together with Proposition IA.2I and Lemma [2.2l . we have the following result. 

Proposition A. 10. Let Dj be as above, m G N and p € M, 1 < p < oo. Assume m > d/p 
when p > 1, and m > d, for p = 1. If u £ W^{'Dj) satisfies u\x = 0, where X is a finite 

subset of Vj, and if the geodesic meshnorm h = hx < "sm^f^"^' then 

\\''J'\\w^{Vj) ^ C'm,k,p,M^'^~^^iV^)~'^~'^''\\''^\\w^{Vj) (A. 11) 

Ml^(v,) < Cm,k,pMR"'"'^'Piy^)~''Mw-(v,)- (A.12) 
Proof. We first apply Proposition lA. 21 with h the Euclidean meshnorm for Exp~,^ X, satisfying 
h < h/Ti < ^Q^i^i < '^^liHSr^ • '^^^'^i equations (|A.ip and (|A.2p hold for u o Exp^^ on 

Dj, with 7£). and d^. replaced by the bounds in (jA.lOp . Applying Lemma 12.21 then gives the 
bounds above. □ 

Theorem A. 11 (Manifold Case). Suppose that Q C M is a bounded, Lipschitz domain 
that satisfies a uniform cone condition, with the cone having radius Rq < rM/3 and angle ip. 
Let k, m, and p be as in Proposition \A.S\. and let X ci^ be a discrete set with mesh norm h 
satisfying 

h<hoRn, ''o:=^^^F{ipf, (A.13) 

where ri,r2 and F(-) are defined in i2. 3\) and II A . 3\) . respectively. If u £ W^{Q) satisfies 
ulx = 0, then we have 

\Mw,Hn) < Crn,k,PMl'i^y^''''^'^''~'"'h"'~'M\wp^in) (A.14) 

and 

hh^^n) < C^,p,M/i"^-''/Pi^(v')-'+''/P-'"^||n||w™(n). (A.15) 
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Proof. We are given h in (jA.lSP to begin with. Thus, we may choose R = < 
and also take the 'Dys to be the domains corresponding to this R. It follows that the 
conditions on h in Proposition lA. 10] hold: consequently, 



2r4 



lkllH/p'=(©j) < Cm,k,p,M^^ ^i'P) "^\\'^\\w^(Vj), (A. 16) 

||n||z.^(^,^) < C™,p,M/^"-"/PF((^)-'^+''^/P-'™||n||H.^^(i,^). (A.17) 

Because of the decomposition = Up^^T^Vj, the bound in (|A.17p immediately implies (jA.lSp . 
Moreover, this decomposition also gives us 



From Definition 12.11 we see that 

^ M\w^{Vj) = ^ j^^Xv,{p)\'^'f\l^p^tJ'{p) < sup iY^Xv.ip) 



I 11^ 



j i=0 



The sum xVj{p) is precisely the number of P^'s that contain p. Since Vj C b(pj,i?), pj 
is itself also in b(p, R). Consequently, the number of ^?/s containing p is bounded above by 
the number of balls b(pj,r), where r = 2RF{ip) ~ h/F(ip), that contain p, and, ultimately 
by the maximum number of b(pj, r)'s that can intersect each h{p, R). By Lemma lA.51 this is 

3(d-l) ^^^^ 

(4a + l)'^e , where a = R/r = ^F{ip)~^. Putting together the two previous inequalities 

then yields 

Taking the p*'' root, lumping constants, and manipulating the result, we obtain ()A.14p . □ 

We remark that the various constants appearing in Theorem lA.lH including /iq, only 
depend on and only the right side of (|A.13p depends on the radius Rq, and that dependence 
is linear. Thus, the dependence on is completely explicit. 

At this point we can extend the Duchon type error estimates for approximation by con- 
ditionally positive definite kernels. To our knowledge, this is the first result of this kind on 
bounded regions in compact Riemannian manifolds. 

To this end, suppose i^T : x 17 M is positive definite (i.e. the matrix 
is positive definite for each H - see Definition 13. Ih and consider that the "native space" 
Mk, the reproducing kernel Hilbert space constructed by taking the space of arbitrary linear 
combinations of completed under the inner product (/,<?) i— )• ^ C) for 

/ = ^ A^K{-,(^) and g = -^C-^(') 0- In this case, it is well known that the kernel interpolant 
Is/ is the optimal interpolant in the sense oi Mr- Namely, ||/=/||a/'a' ^ I|s||a/'x for s G Mr 
with s|h = = . 
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Corollary A. 12. Let m > d/2, and let K be a positive definite kernel on M for which Mk is 
continuously embedded in W2^{Wl). Let C M satisfy a uniform cone condition with radius 
Rci < I'm/3 and angle ip. For H C having meshnorm h < hoR^i and for f G Mr, 

We note that this result holds for a much larger class of kernels than considered in the 
previous sections (i.e., defined by Definition 13. 2p . In particular, there are numerous examples 
of compactly supported kernels on M*^ and S*^ having native spaces that are Sobolev spaces, but 
which do not invert an elliptic differential operator. However, his type of error estimate should 
be compared to those in Corollary 15.91 and Corollary 15.101 - observe that the condition on the 
target function is quite restrictive (it needs to be in Mr and there is a basic disagreement 
between the approximation order m — d/2 and the smoothness assumption). 

Proof. Apply Theorem I A. 11 1 to u = f — I^, we see that 

ll/-WllL(n) < 
< 

< 
< 

The next to last inequality is the embedding Mr C W2"'{M.), while the last inequality is the 
Pythagorean theorem for orthogonal projectors ||/ — IsfWxj^ H-^^/Ha/'k ~ II/Ha^k ' 

There are several domains that are important for us, and that we will discuss below. 
We begin with the manifold itself. In that case, we may take Rn = rM/3. The angle <p 
may be set equal to vr/2, because every such cone is contained in M. This means that 
F{ip) = F(vr/2) = 1/8. 

Corollary A. 13 (Full Manifold). Suppose that M is compact. Let k, m, and p he as in 

p2 

Theorem \A . 1 1[ Then, with ho := gs^lpij , there is a constant Cm,k,d,M such that i/ A C M has 
mesh norm h < /iorM/3 and if u G W^{M) satisfies u\x = 0, then 

\\'^\\w^{M) ^ C'm,fc,p,M^'" '^ll'"lll4^p™(M)- (A. 18) 

The domains that we now turn to are balls, annuli, and complements of balls. In all of the 
cases discussed below, the domains Vt satisfy the ball property described in Corollary IA.8I at 
each point q G dVl. Consequently, we may take (/? = arcsin(2^), and so F{(p) = l/(8r2 + 4). 
It thus follows that in all such cases 

"'^ = 128m4k + l)^ - ^^-^^^ 



{crn,Mh'^-'^'Fi^r''^y 11/ - m'w-m 

C (c^,m/i'"-''/''F(<^)-2-)' 11/ - L^fWl^ 

CKF{^r'^h^^^-WNK- 
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and all of the factors in Theorem lA.llI depend only on parameters from the manifold itself, 
as well as p,k,m, but not at all on the center and the radius of the ball/annulus/punctured 
ball. 

We now turn to balls. A ball of any size may be treated; however, if the the radius is larger 
that rM it may intersect itself, giving rise to corners with angles that have to be dealt with 
on a case by case basis. This phenomenon is easy to see in the case of the torus embedded in 
M^, where a sufficiently large ball begins to wrap back on itself. When the radius of the ball 
is less than rM, this wrapping doesn't happen. With this assumption, we have the following 
result: 

Corollary A. 14 (Zeros estimate on balls). Assume m > d/2. Suppose that r < tm- If 
u G ty™(M) vanishes on X C b(p, r), where h < h^r/l, where /iq is given by liA.19\) . we have 

II II 7 TTh A; I I 

Proof. Obviously every point q £ dh{p,r) satisfies the conditions in CoroUarv IA.8I A direct 
application of Theorem lA.llI then completes the proof. □ 



Corollary A. 15 (Holder estimate on balls). // m is greater than d/2 + e, and the 

conditions of Corollary \A.14 'iold (in particular, r is less than tm, cind h < h^r), and if 
u G Ty™(b(p, r)) satisfies u\x = 0, then for every z G b(p, r), 

\u{p) - u{z)\ < C7r™-^-'^/2^ist(p,z)^||n||pvm(b(p,^)), 

where C is a constant depending only on m, M and e. 

Proof. This follows because the Sobolev embedding theorem, in conjunction with Lemma 
ensures that w = w o Expp € C^(i?(0, rRf))- Thus for z = Expp(x) € b(p, tm), 

\w{p) — w{z)\ \w{Q) — w{x)\ 



< 



w\c-{BiO,rM)) < C\\w\\w^(B{0,tm))- 



dist(p, z)"^ Ixl*^ 
For a general r < rM, set ui^^x) = u{x). Then 




1/2 



The result follows by applying Lemma 12.21 in conjunction with Corollarv IA.14I □ 

A similar argument to the proof of Corollary IA.14I given above yields these results for annuli 
and complements of balls. In the following two lemmas, we are concerned with the case p = 2. 
Consequently, we suppress dependence on these parameters by expressing the constant from 
the zeros lemma for such domains simply as A. In other words. 



v{3/2)d+2m 

iiia,^ '-'r?i,fe,::i,M V"-^ ^ ^ 
k=0...m-l 

which depends only on m and 



A:= max C^,fc,2,M (SPa + 4)^'/'^'*+'" , (A.20) 

k=0...m-l . ' . ^ ' 
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Corollary A. 16 (Zeros lemma on annuli). Assume m > d/2, Let a = h{p, r) \ h{p, r — t), 
where < i < r < im, and let /iq he given by iA.19\) . If u £ W^(a) vanishes on X C a., 
where h < /iq min(i/2, rM/3), we have 

ll^llw|(a) ^ ^h"^~'^\u\w^{a)- 

Proof. At each point q of the boundary of a an open ball of radius t/2 can be placed inside a 
with a boundary that passes through q. The result follows from Corollary IA.8I and Theorem 

ixm □ 

Corollary A. 17 (Zeros lemma on complements of balls). If r < rM/3 and if u £ 
VF™(M) vanishes on X with h = h{X,h{p,rY) < /iorM/3, then 

Proof. By placing its center, q, a distance of r + 2rM/3 away from p, the ball h{q, 2rM/3) can be 
placed in h{p, rY- It follows that the set satisfies a cone condition with radius R = rM/3. □ 
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